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Abstract 

We study the dynamics of the populations of a model molecule endowed with two sets of 
rotational levels of different parity, whose ground levels are energy degenerate and coupled 
by a constant interaction. The relaxation rate from one set of levels to the other one has an 
interesting dependence on the average collision frequency of the molecules in the gas. This is 
interpreted as a quantum Zeno effect due to the decoherence effects provoked by the molecular 
collisions. 

Introduction 

The quantum Zeno effect is usually formulated as the hindrance of the evolution of a quantum 
system due to frequent measurements performed by a classical apparatus-^ and is formalized 
according to von Neumann projection rule. 3 The literature of the last few years on this topic is vast 
and contemplates a variety of physical phenomena, ranging from oscillating (few level) systems 4 
and alternative proposals- to bona fide unstable systems, 6 where the so-called "inverse" Zeno 
effect can take place. 

The ideas and concepts at the basis of the quantum Zeno effect (QZE) were also successfully 
extended to continuous measurement processes by different authors and in different contexts 7 and 
led to a remarkable explanation of the stability of chiral molecules. 8 This was a fertile idea, in 
that it explained the behavior of a variety of physical systems in terms of a similar underlying 
mechanism. 

The QZE is, however, a much more general phenomenon, that takes place when a quantum 
system is strongly coupled to another system- or when it undergoes a rapid dephasing process. 
Such a rapid loss of phase coherence ("decoherence") of the quantum mechanical wave function 
(for instance as a result of frequent interactions with the environment) is basically equivalent to 
a continuous measurement process (the main difference being that the state of the system is not 
necessarily explicitly recorded by a pointer). 
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The quantum Zeno effect is always ultimately ascribable to the short-time features of the dy- 
namical evolution law:— it is only the study of this dynamical problem that determines the range 
in which a frequent disturbance or interaction will yield a QZE. The very definition of "frequent" 
is a delicate problem, that depends on the features of the interaction Hamiltonian. Moreover, one 
should also notice that the quantum system is not necessarily frozen in its initial state, 11 but rather 
undergoes a "quantum Zeno dynamics", possibly evolving away from its initial state. 13 The study 
of such an evolution in the "quantum Zeno subspace" 12 is in itself an interesting problem, whose 
mathematical and physical aspects, as well as the possible applications to chemistry and physical 
chemistry, are not completely clear and require further study and elucidation. 

After the seminal experiment by Itano and collaborators, 4 the QZE has been experimentally 
verified in a variety of different situations, on experiments involving photon polarization,— nuclear 
spin isomers, 15 individual ions,-^"— optical pumping,— NMR,— Bose-Einstein condensates,— 
the photon number of the electromagnetic field in a cavity, 23 and new experiments are in prepara- 
tion with neutron spin.- 2 ^ 2 ^ 

We focus here on the interesting example of QZE proposed in:— the nuclear spin depolarization 
mechanisms in 13 CH3F, due to magnetic dipole interactions and collisions among the molecules 
in the gas, was experimentally investigated and interpreted as a QZE. In a few words, the 13 CH3F 
molecule has two kinds of angular momentum states, according to the value of the total spin of 
the three protons (H nuclei): 1 = 3/2 (ortho) and 1=1/2 (para). Transitions between states with 
different parity are (electric dipole) forbidden, so that spin flip occurs via a weak coupling between 
two levels of different spin parity (this is most effective when there is an accidental degeneracy 
between the levels, achievable, for example, via a Stark effect—). One observes a significant 
dependence of the spin relaxation on the gas pressure and interprets this as a QZE provoked by 
the dephasing due to molecular collisions. Nuclear spin conversion in polyatomic molecules is 
reviewed in.— 

The aim of this article is to study the occurrence of the QZE in the general framework of 
collision-inhibited Rabi-like oscillations between two sets of rotational levels. We shall study 
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the evolution of the level populations in a model multilevel molecule endowed with two sets of 
rotational levels of different parity. In particular, we shall concentrate on the interesting effects 
that arise as a consequence of the interactions (collisions) with the other molecules in the gas. 
The model we shall adopt will be studied both numerically and analytically, and the results will 
be compared. One of the main objectives of our investigation will be the analysis of apparently 
different phenomena in terms of a Zeno dynamics. 

We shall introduce the system in Section„and the Zeno problem in the present context in Section 
„ In Sections □-□ we study the problem from an analytic point of view, by deriving and approximately 
solving a master equation. In Section „ the analytical result are compared to an accurate numerical 
simulation. We conclude in Section „ with a few remarks. 



The system 

Our model molecule has two subsets of rotational levels (to be called left (L) and right (R) levels 
in the following) of different parity, whose ground levels are energy degenerate and coupled by 
a constant interaction. (The choice of the ground levels is motivated by simplicity: one could 
choose any other couple of energy-degenerate levels in the L-R subspaces.) The molecules undergo 
collisions with other identical molecules in the gas and we assume that these collisions couple 
the rotational energy levels but do not violate spin parity conservation. We shall focus on the 
dependence of the relaxation rate on the average collision time or, equivalently, on gas pressure: 
a QZE takes place if the transition between the left and right subspaces is inhibited when the 
collisions become more frequent (i.e., the gas pressure increases). 



A sketch of the system is shown in Figure 1 The total Hilbert spaces of each molecule is 



made up of two subspaces Jtft (left) and J$? R (right), with Nl and levels respectively. Collisions 
cannot provoke L <H> R transitions, so that no transitions are possible between the two subspaces, 
except through their ground states. However, collisions with other particles in the gas provoke 
transitions within each subspace. 
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The Hamiltonian is 

H = H f + H con (t)=H + H 1 +H coll (t) 1 (1) 
where Hf =Hq + H\ is the free Hamiltonian and 

N L N r 

Hq = E n L \ n L){ n i\+ Y E nR \n R )(n R \, (2) 

n L =l n R =l 

H x = Mi(|l L )(l*| + |l*)<l L |), (3) 

#coll = T/)V, (4) 

i 

V = (XlVl + ccrVr, (5) 

Ns-l N s -i 

V s = Y, V »s = L \n s ){n s + l\ + \n s + l){n s \, (6) 

n s =l n s =l 

with s = L,R. The energy levels \n s ) have energies E„ s (s = L,i?) and H\ provokes L <H> R transi- 
tions between the two ground states, with (Rabi) frequency Q.. £1 is small (in a sense to be made 
precise later), for such a transition is electric-dipole forbidden. H co \\ accounts for the effect of col- 
lisions with the gas (environment): the collisions are distributed according to the Poisson statistics, 
appropriate for gas phase with short-range binary collisions, so that they occur at times 

T i+i = t ; + 5t ;' ( 7 ) 
where 5t/s are independent random variables with distribution 

p(8xj) = -e,xp(-8xj/r) (8) 

and (common) average x. The coupling constants OLl r are in general different from each other 
and measure the "effectiveness" of a collision. For the sake of simplicity we assume that colli- 
sions provoke transitions only between adjacent levels [V s in © involves only "nearest neighbors" 
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couplings]. We will assume, for concreteness, that the energy levels are rotational, so that 



E ns = h(D s n s {n s + 1) (s = L,R) (9) 
and | lx) and \1r) are the only resonant pair of states: 

E\ L = E\ R , E mL ^E llR for m L ,n R >l. (10) 



See Figure 1 The Hilbert spaces Jft and ,Jrf?R are finite dimensional, with dimensions Nl and 
Nr, respectively. This is because, in general, the number of accessible rotational levels is limited 
to a few tens, since for sufficiently high energies molecules tend to dissociate. This could be 
accounted for by introducing two "absorbing" levels \Nl+i), \Nr+\).— However, in our analysis, 
we will explore a time region in which the introduction of absorbing levels is not necessary (in 
other words, the times involved will not be long enough to display "border effects"). 

Zeno effect 

Before we start our theoretical and numerical analysis it is convenient to focus on the physics of 
the model introduced in the preceding section and to clarify in which sense we expect a Zeno effect 
to take place. We start from a simple numerical experiment and calculate the time evolution of the 
populations by the Monte Carlo method described in.— 

Consider a uniform gas of identical molecules, having the internal structure described in the 
preceding section. A single molecule freely wanders in a total volume and undergoes random colli- 
sions. By neglecting the spatial component of the wave function, each molecule can be represented 
by an (Nl + AO;)-dimensional state vector | that describes its internal state. 30 This physical 
situation is well schematized by the model described in Sec. a During the free flight the evolution 
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is governed by the free hamiltonian 



| !//■(*)) =exp 



Mo)). 



(ii) 



Since the molecules are immersed in a bath, the collisions are distributed in time according to 
the Poisson statistics ©-©, with average collision frequency (per particle) T . Once a collision 
occurs, a collision time is sampled according to: 



y being a random number uniformly distributed in [0, 1 [. The collisions are modelled as instanta- 
neous events and act on the left/right subspaces independently. As a result of a collision, the state 
becomes 



The matrix exp(— ?X S 0C s V s ) is evaluated numerically. It is assumed to be independent of the internal 
state of the colliding partners and of their kinetic energy. 

We also stress that since our aim is to investigate the occurrence of a QZE within the proposed 
level structure, we are not interested in the dissociation of highly excited molecules. To this end, 
we must restrict our attention to times such that the molecules do not "see" the upper limit of 
the rotational levels, so that "border" effects do not play any significant role. In this way the 
dissociation of highly excited molecules can be safely neglected. 

The afore-mentioned qualitative features of our analysis will be carefully scrutinized and made 
precise in the following sections. We now take them for granted and give a few preliminary results 
in order to get a feeling for the physics at the basis of the Zeno effect. 

We set Nl = Nr = 40 energy levels, with energies given by ©, where n s = 1, . . . ,40, (Oi — 
1.3 • 10 10 s _1 and (Or = 9.7 ■ 10 9 s~ ! . We always compute the average over an ensemble of 5 • 10 3 
particles. All particles are initially in the \ l L ) state and we study the temporal behavior of the 



8x 



Tlog(y), 



(12) 




(13) 
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relative population in the left subspace 



(14) 



I'L 



p„ L being the occupation probability of state \ni). 



The results of our numerical integration are shown in Figure 2 Figure 3 In Figure 2, ai — 0.2 



1,00 



oc L =0.2 , a R =0 
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Figure 2: Temporal evolution of Pi. The collision frequency T _1 is varied between 5007^" and 
15007^ (T R = 2n/£l). We set a L = 0.2, a R = 0, so that, in practice, N L = 40 left energy levels 
are coupled to only Nr — 1 right level. The survival probability in the left subspace increases as 
the collision frequency is increased: frequent collisions hinder transition to the right subspace, a 
manifestation of a ("classically intuitive") Zeno effect. 



and (Xr = 0, so that collisions do not provoke transitions among the right states (or, equivalently, 
the right subspace consists only of state |1^)). It is apparent that when the collision frequency 
T _1 is increased between 5007^ 1 and lSOOJ^ 1 (Tr = 2tz/Q. being the Rabi period) the survival 
probability in the left subspace increases. If the collisions are viewed as a dephasing process 



9 
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Figure 3: Temporal evolution of Pi. The collision frequency % 1 is varied between 500T R 1 and 
lSOOT^T 1 . Unlike in the previous figure, we set at = 0, CCr = 0.2, so that in practice, Nl = 1 
left level is coupled to = 40 right levels. Again, the survival probability in the left subspace 
increases as the collision frequency is increased: frequent collisions hinder transition to the right 
subspace, a manifestation of a ("classically counterintuitive") Zeno effect. 
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(effectively yielding a "measurement" of the occupation probabilities of the left states), this can 
be viewed as a Zeno effect. This is in agreement with our "classical" intuition: since the system 
is initially in the left subspace and collisions remove population density from the ground state 
of this subspace (the only level coupled to the right subspace), it is intuitively clear that, 
by increasing the collision frequency, transitions to the right subspace are hindered. This is a 
"classically intuitive" version of the Zeno effect. 



The situation depicted in Figure 3 is different: here (Xl = and a# = 0.2, so that now collisions 
do not provoke transitions among the left states, or equivalently the left subspace consists only of 
state (which is also the only state coupled to the right subspace). Once again, when the 
collision frequency T _1 is increased in the same range as before, the survival probability in the left 
subspace increases. We stress that the collisions are effective in hindering the transition from a 
single level towards a subspace that is initially empty. In other words, now the collisions act only 
on the right subspace, where virtually no particles are present. Once again, this can be viewed as a 
Zeno effect; however, it is somewhat less intuitive than the previous one (and maybe a bit puzzling 
for our "classical" intuition). This a "classically counterintuitive" version of the Zeno effect. 

After having rapidly analyzed these two simple situations, we are ready to tackle the more 
general case of Nt left levels coupled to right ones. This will be done in the following. 

Master equation 
The general case 

We start our analysis by deriving a master equation for the density matrix of the molecule. Write 
© as 

#cqu(0 = M0V, (15) 



where 

, \ d 
nit) = - 

dt 



M(0 = ^ = E*(f-*y) 06) 



n 



is the derivative of a Poisson process N(t) with mean time t:— 

Ptob{JV(0 = n} =P(n,t) = e~ tlx — l (-Y . (17) 

One gets 



so that the process 



(dN(t)) = d l l(dN(t)- d A 2 \ = d l (18) 



dW(t) = r] (t)dt = n(t)dt - — = dN(t) - — (19) 

T X 



has a vanishing mean and a linear variance in dt 



dt 

(dW(t)) = 0, (dW(t) 2 ) = — . (20) 

x 



In terms of the white noise Tj (t) these equations read 



= 0, (t?(Ol(0> = ^(*-0- (2D 



The collision Hamiltonian can then be rewritten in terms of a constant part and a white noise 

H con (t) = -V + hr 1 (t)V, (22) 

T 



whence the total Hamiltonian (OQ) reads 



h 

H = H + hr)(t)V, H = H +H 1 + -V. (23) 

x 



The Schrodinger equation (in Ito form) is 

i - 



\}/(t + dt)}= ( l--Hdt-^-V 2 dt ) \ Y(t))-iVdW\\ir(t)} (24) 
n 2t 
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and the average density matrix [(•■■) is introduced in Eq. (fT8l)1. 



p(t) = (\Y(t)){xir(t)\), (25) 



follows a master equation in the Kossakowski-Lindblad form 



32 



= - l j i \fl,p]-^{V\p}+ l -VpV. (26) 



By using Eqs. d2~3l . © and © we get 



f - 4^p]4^p]4[v,p]4{v 2 ,p} + ^pv 



\ [H ,p] -^[H h p]+ ^ ("^,p] - ^{V^P} + ^V sP V s 



+^-(V L pV R + V R pV L ). (27) 



This equation for the average density matrix (1231) is exact, but complicated. However, it can be 
greatly simplified under some reasonable hypotheses. 

Reduced master equation 

We assume that all level pairs E mR and E nL are sufficiently far from resonance, namely 

AE , 

— »Q, T'\ (28) 
n 

where AE is the smallest energy difference between states |m^) and |nj,), with mR.ni > 1. This 
requirement will be discussed in more detail in Section „ At this stage we only observe that, 
typically, AE/H ~ 10 -9 s, while Q. ~ 1kHz and % ~ 1/is, so the above condition appears very 
reasonable. 
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It is then possible to show that in (1271 ) the dynamics of the populations p Ms = p m ,m s (s = 
L,R) plus the coherence term P\ L \ R completely decouples from the dynamics of the coherence 



frequency is present [essentially because (m s \[Ho, p]\m s ) = in (l27l) l and, under hypothesis (1281) . 
the contribution of all the other fast terms is averaged to zero over the long timescales x and 
and the dynamics of the slow and fast terms completely decouples. In conclusion, only the "slow" 
dynamics is relevant over the large timescales x and Qr x . 

The above argument has a general rigorous justification— in terms of an adiabatic theorem and 
is elucidated in Appendix n for the model studied in this article. One shows that the part of the 
master equation ([27l) pertaining to the populations becomes 



terms p msn , (s,s' = L,R and m s ,n s i ^ 1). This is because, roughly speaking, no "diagonal" fast 



dp 



[H h p)-^[v,p)-^{V^p} + ^VpV, 



(29) 



dt 



where the reduced operator A, defined by 



A = QAQ+ £ £ P ms AP, 



(30) 



s=L,Rm s =2 



involves only matrix elements belonging to the eigenspaces of Hq,— 



Q = Pi l +Pi 



R 



and P, 



m s — 




(3D 



[remember condition (fl0l) l and is diagonal with respect to Hq 



[#o,A]=0. 



(32) 
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In particular, from Eq. © 



Hi = QH 1 Q = H l =hao h 

01 = + (33) 



and from Eqs. © and © 



and 



y = o, v 2 = £ a 2 V 2 = £ a 2 "£ V 2 (34) 

s=L,R s=L.R fj!,= l 



ypy = £ «2y,py = £ « 2 £ y m . s py m „ (35) 

s=L,R s=L,R m s =i 



so that Eq. (1291) reads 



^ = -in[a 1} />] - £ £ £ [y ms , [y„,,p]]. (36) 

at s=L,R ZT m*=l 

This is the master equation we will study in detail. The only assumption made in its derivation is 

The reduced density matrix p is given by Eq. (1301) and involves only the level populations 
Pn s = pn s n s (s = L,R) and the two coherence terms pi L) i s and pi s ,i L , all other matrix elements 
being zero. Thus, it describes two classical Markov chains (1l, ■ ■ • ,Nl) and (1^, . . . ,Nr), whose 
transition rates are proportional to D L = 05 2 /t and Dr = a\/x respectively, linked by quantum 
Rabi oscillations between | \i) and \ \r), whose period is Tr = 2n/Q.. 

By setting D = (Di + Dr)/2, i.e. a 2 = (a| + a|) /2, the ratio between the two timescales 

_ oc 2 2na 2 
x = DT R = —T R = — — 37 

is an important parameter, that describes different dynamical regimes. Larger values of x cor- 
respond to more frequent collisions (within a Rabi period) and consequently to a more evident 
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manifestation of the QZE. 



Stochastic dynamics in decoupled subspaces 



Let us first study the subdynamics of each subspace separately. To this end, set £1 = in Eq. 
(l36l) : the time evolution is governed only by the collision dynamics, the right and left subspaces 
decouple and their subdynamics can be studied separately. 
In terms of (s = L, R) 



P S = (PlAvjPn^,...) = (Ph,---,Pn s ,---) 



(38) 



Eq. (1361) reduces to 



where W s is the stochastic matrix 



dp s 
dt 



D s W s p s 



(39) 



and 



/-ll 
1-210 
1-21 
1-2 

V : : ; ; 



«7 
T 



(s = L,R) 



(40) 



(41) 



Note that W s is a real symmetric matrix with real eigenvalues and a complete set of eigenvectors. 
The resulting dynamics is diffusive. Indeed, Eq. (l39l) explicitly reads 



p ls = D s (-p ls +p 2s ), (Q = 0) 

Pn s = D s {pn s -\-2p ns +p ns +\), (n s >2) 



(42) 
(43) 
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which is nothing but a diffusion equation (dropping the suffix s) 

d t p n (t)=DAp n (t), 



(44) 



where 



A = I(V+V~+V - V + ), (45) 



and 

V + p n = Pn+1- Pn, V~Pn = Pn~ Pn~\- (46) 

The boundary condition V~pi = [see Eq. (1421) 1 is imposed by introducing a supplementary state 
n = 0, whose probability satisfies po = p\ for every t. The evolution of the population is made 
up of two terms (both expressed in terms of the fundamental solution of the heat equation): each 
"site" (level) gets a direct and a "reflected" contribution from the boundary n= 1: 

p n (t) = q n (Dt)+ qi - n (Dt), (47) 

where q„(t) are the probabilities of the continuous-time symmetrical random walk engendered by 
the equations 33 

q n = q n -i -2q n + q n +l, -°°<n<oo 5 (48) 
whose solution starting at n = 1 for t = [i.e. p n (0) = 8 n i] reads 

q n (t) = e- 2t I [n _ ll (2t), (49) 

I n (t) being the modified Bessel function.— As is well known, for t — > °° and n — > °o with n 2 /t =const, 
this yields 

g„(0~— Lexp f- ^T/H - ( 5 °) 



/4m V 4? 
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Note that the boundary condition is essential in assuring probability conservation, 

-\-oo -\-oo 

!>»(*)= E q n (Dt) = i, (5i) 

n=l n=—°° 

for any t. The above equations are of general validity. In particular, 

Pl (t) = e- 2Dt [I (2Dt)+h(2Dt)} , (52) 

(and pi(0) = 1). It is also possible, by using the solution (|47T) and (I49)) . to evaluate the mean and 
second moment 



H(t) = £n/>„(f), (53) 

«=1 

CO 

a 2 (t) = (54) 

Indeed, by using (1481) . one can obtain explicit differential equations involving these quantities, 
valid for any t, 

fl(t) = D Pl (t)=De- 2Dt [I (2Dt)+h(2Dt)], 

o 2 (t) = 2D+D Pl (t), (55) 

whose integration gives 

H(t) = l -+ l -e- 2Dt [{\+ADt)I Q {2Dt)+4Dth{2Dt)}, 

(56) 

o 2 {t) = 2Dt + ti{t). (57) 
Let us also give, for completeness, the expression of /i and o for times Dt 3> 1 [but always 
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t < T d , see Eq. <[60]> below]. From Eqs. (I56l)-(I57T). 



[AD 




(58) 



In order to compare these results with those of the following sections, consider that 



t 

= x — , 
Tr 



(59) 



where T R is the Rabi period and x = x s — ci 2 Tr/x (s = 1,2) is essentially the scaling parameter 
introduced in dTTl) . 

In reality, as we already emphasized, the number N (Nl or A^) of accessible rotational levels 
is in fact finite, because the molecule dissociates after absorbing a sufficient amount of energy. 
In order to account for this process one can add an (N + l)th absorbing level (in each subspace). 
However, since we are interested in phenomena, such as the QZE, that can be brought to light 
within timescales shorter than the dissociation time, the introduction of an absorbing level is an 
unnecessary complication that can be easily avoided by restricting our attention to the relevant 
timescales. Let us therefore estimate the time scale at which dissociation occurs. If a molecule 
dissociates when it reaches level N + 1, namely if only N levels take part in the diffusion process, 
the "dissociation" time reads 



This is the time needed by the system, that starts in the ground level, in order to reach the upper- 
most level via the diffusive propagation engendered by the collision. This rough estimate of the 
"dissociation" time can be improved: a better analysis yields — N 2 x/n 2 a 2 , which is roughly of 
the same order of magnitude. 

In our analysis we will assume N = 40. Within the numerical range of the parameters a and x 
to be used in our simulation, the dissociation time 7^ varies between 2 and 12 Rabi periods. In the 



N 2 x N 2 



(60) 



Td - — 



D a 2 x 
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following we will always remain well below this threshold. 

Zeno effect in coupled subspaces 

We have seen in Section n that when the L and R subspaces are coupled, namely when Q. ^ 0, a QZE 
can be obtained by increasing the collision frequency. Indeed, as we will show, by increasing the 
collision frequency, the probability of remaining in the initial subspace decays more slowly. We 
also commented on the possibility of studying the Zeno dynamics in two different situations, one 
classically more intuitive and the other one less intuitive. These different names reflect the fact that 
the former case can be understood (at least qualitatively) by means of a classical Markov process, 
while the latter cannot. Both Zeno effects are contained in the master equation derived in Section 
n and are a consequence of the features of the collisions with the other molecules constituting the 
environment, or in other words, of the coupling constants of the interaction Hamiltonian H co \\. The 
resulting dynamics will be numerically investigated in full generality in Section „ However, before 
we show the results of the numerical simulation, let us discuss the main qualitative features of the 
dynamics without solving the complete master equation. This will be done in the present section 
with the help of some working hypotheses and will help us clarify some additional features of the 
Zeno effects. 

When Q ^0 the two subspaces Ml and Mr are coupled through their ground states. The 
evolution is described by (l43l) . supplemented by the following three equations 

pi L = D L (-p lL +p 2L ) + ap c (t), (61) 
Pl R = D R (-p lR +p 2R )-a P c (t), (62) 
fit) = -^±^p^t)-2a(p lL (t)-p lR (t)), (63) 

where p c is the coherence term between states 1 1^) and \1r) 

p° = -2Impi L i fl = i(pi L i R -p\ R \ L ) , (64) 
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responsible for the coupling between the two subspaces. The total probabilities of being in the left 
and right subspaces read 

ft(0 = Ep» i (0, fle(0=£^(O s (65) 

nL m 

respectively. The derivatives of these quantities are easily seen to be simply related to the coherence 
term: 

P L (t) = a P c (t), P R (t) = -£lp c (t). (66) 

Notice that Pl + Pr = (conservation of particles number). Let our particles start in the left sub- 
space at time t = 0. Therefore the quantity of interest is P^. One can obtain the evolution equation 
for by eliminating p c by means of ( |63l ) 

P L + DP L + 2Q. 2 (p lL -p lR )=0, (67) 

where we set D = (Dl + Dr)/2. This equation shows that the dynamics of Pi is governed only by 
the population difference between the ground states, irrespectively of the population of the higher 
levels. This introduces an interesting picture of the dynamics, in which the Rabi oscillations act 
as a "source" for the probability. The source drains particles from the left to the right subspace if 
p\ L > p\ R and vice versa if p\ L < p\ R . 

Our initial condition will always be p\ L {0) = 1: initial population in the ground state of the 
left subspace. For t <C \/D <^Tr = 2tz/Q. (which is always true for our choice of parameters), we 
can set p\ L {t) = 1 + O(Dt) and p\ R (t) = O(Dt) and a power-series solution of (1671) . with initial 
conditions P L (0) = 1, P L (0) = 0, yields 

P L {t) = 1 - QV + o(Q. 2 t 2 ) , (68) 

which shows that the quadratic region of the Rabi oscillation is not perturbed by the collisional 
dynamics (namely, does not depend on a), even thought it extends up to times shorter than 1/D 
Tr. This result was to be expected 9 and is well observed in our numerical experiments, also for 
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very high collision frequencies. 

Equation (1671) is exact, but it is not a closed equation for the total probability P L . One needs the 
populations of the ground states in order to obtain Pi. We will therefore introduce an ansatz for 
the functional form of the populations of the ground states, valid for large values of the parameter 
x defined in Eq. (1371) . which will enable us to get a closed equation for Pi. In addition we will 
also gain a deeper understanding of the Zeno phenomenon for this system. The ansatz consists in 
substituting for p\ L/R the solution (l52l) . obtained for the decoupled subspaces (Q. = 0), normalized 
to Pi/r. This "adiabatic" (Born-Oppenheimer-like) approximation relies upon the assumption that 
the time scale of the internal collisional dynamics is much faster than the Rabi one (1/D «C T R ), 
so that particles are drained from the ground level and redistributed according to the uncoupled 
dynamics. The Rabi coupling simply accounts for the varying number of particles present in each 
subspace. This ansatz is translated into the equations 

Pl L = P L (t)f L (t), 

PU = P R (t)f R (t) = (l-P L )f R (t) 7 (69) 

where fi/R(t) are the population probabilities of the ground states given by the uncoupled dynam- 
ics 452} 

k/R = e- 2DL ' Rt [k(2D L/R t) + h (2D L/R t)] . (70) 
Substituting in (1671) we obtain 

Pi + DP L + 2a 2 (P L (fi + /*)-/*) = 0, (71) 

which is the equation of motion of a unit-mass, forced pendulum with varying frequency. The 
initial conditions are Pl(0) = 1, Pl(0) = 0. It is easy to prove that if f R / (fjL + fit) tends to a well- 
defined limit and its first and second derivatives vanish when t — > °o, there is a stable fixed point at 
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t — 



oo 







/*(') 



(72) 



fR(t)+f L (t) 



and any solution will eventually reach this point. This feature of the population of the left subspace 
is always observed in the numerical solutions. An asymptotic analysis of the Bessel functions, 
performed with l/D <C t <C 7d shows that all these requirements are satisfied and an equilibrium 
distribution exists and is given by 



Let us see now how the Zeno effect emerges in this picture in three different cases. 

Case ccl = ccr = ccj^O 

The first case-study is obtained by setting a = OLl = (Xr ^ 0, so that the (collision dynamics in the) 
two subspaces are identical and f L (t) = _//?(?) = f(t). We change the time variable from t to the 
dimensionless t /Tr — 2nt/Cl and set x = DT R = <x 2 Tr/t = 2na 2 x/Q., obtaining (the dot denotes 
now differentiation with respect to t / Tr) 



where f(t) is given by d70l) . Since D L = D R , according to (1731) . Pl will eventually tend to P[ = 1/2. 
However, we will see that the typical time scale r re i ax of this relaxation process will increase with 
x and this will be interpreted as a QZE. 

The proposed analogy with a classical damped harmonic oscillator suggests that when x 3> 1 
we get Pl ~ 1 /x a and Pi ~ 1 /x 2a , with a > 1. Indeed, we will see that the solution satisfies this 
hypothesis with a = 3/2, so that the first term (Pl ~ x~ 2a ) is negligible with respect to the second 
(xPl ~ x l ~ a ) and the third one (both / and Pl are of order 1) and hence can be dropped from (l74l) . 

'Actually one always (implicitly) assumes t <C 7d- If one looks at longer times t > T^, this equilibrium state appears 
as a metastable state, which decays into the true equilibrium state. 




(73) 



1 + ^JDr/Dl ocl + ccr' 



PL+xP L + S7t 2 f(t)(2P L -l) =0, 



(74) 
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Thus we are left with a first-order, separable differential equation whose solution is [here Pl(0) = 1 
butP L (O) = 0(l/x)] 



Pdt) 



1 1 

2 + 2 eXP 



Sn 2 



(l-e- 2xt (l+4xt)I (2xt) 



-e- 2xt 4xth(2xt] 



(75) 



For if > 1 we obtain a stretched exponential 



^ 1 1 
Pdt) ~ - + -exp 



Sn 2 



■32n 3 / 2 



t \ 1/2" 
^3 



(76) 



from which one can define a relaxation time as the only characteristic time present in the exponen- 
tial (restoring natural time units): 



Trelax * Tr. 



(77) 



The Zeno effect consists in the fact that by increasing x (more frequent collisions) the correspond- 
ing curves of Pl tend to zero more slowly. These predictions are in qualitative and quantitative 
agreement with the numerical simulations of the next section. 

In order to get a rough preliminary idea of the issues discussed in this section, look for instance 

at 



Figure 4 and Figure 5 , where the numerical results (to be described in greater details in the 



following) are compared to Eqs. (I76l)-(r771). The probability (|75T)-(|761) is correct up to a precision 



of 10%, showing that the ansatz (1691) yields sensible results. Notice that x = 48 in Figure 4 



so 



that the solution (1761) . which is supposed to be valid for xt ^> 1, must yields accurate results for 
t/x 3 > 1(T 6 , as one indeed observes. A numerical fit for the exponent in the stretched-exponential 
yields f 0-3 rather than t 1 / 2 , confirming the general functional dependence. The very fact that the 
global relaxation law is of the stretched-exponential type suggests that the dynamics is highly 



nontrivial, but we will not elaborate on this here. Finally, as can be seen from Figure 5 the scaling 
law (1771) is very well verified. 
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x=48 




Case oll ^ ocr = 



Let us briefly reconsider the first case analyzed in Sec. „ Figure 2 Here the left subspace is affected 
by collisions while the right one is not. Although this is not a realistic situation, it is interesting and 
instructive to look at it. We shall show that also in this case, as the collision strength is increased, 
the system tends to spend more time in the initial (left) subspace. 
If a R = and a L ^ then f L = f and f R = 1 and Eq. (|7TT) reads 

P L + xP L + Sn 2 [P L (l+f(t))-l}=0, (78) 

where x = 7^/21 = Tta^/ Q.T. By means of the same approximations of the preceding subsection 
we obtain, for x 3> 2y/2n, 

%n 2 t i 2V2i N 



1 2^2*-^+ 
P L {t) ~ 1 e 



S7t 2 t 2V27T . 

x* W + , (79) 



where 4>(z) is the error function of imaginary argument— 



<&(z) = -^= \ Z dx/ . (80) 



y/n Jo 

Here the definition of a relaxation time is not easy (no simple scaling law exists). However, both 
in this solution and in the numerical data, Pl has a single minimum P[, which is an increasing 
function of x: this can be regarded as a manifestation of a (classically intuitive) Zeno effect, as 
explained in Sec. „ From (l79l) the value of the minimum is 



2 7 

Pl = l-— (81) 



x 



and is an increasing function of xo This law is well confirmed by the numerical results shown in 



2 Consider f(y) — e y2 <I>(y). Then the numerical value 2.7 in (fSTT > is given by 2y/2nf* , where /* = 0.621 is the 
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Figure 2 



Beyond the minimum Pi tends to 1 with a power-law 



P L{t)-l~\ — ■ 
" TtXt 



(82) 



This is again a Zeno effect: by increasing the collision rate x the survival probability increases. 



Case (Xr ^ (Xl = 



This is the second case analyzed in Sec. „ |Figure 3| If 0^ = 0, (Xr ^ 0, then fa =f,fi, = and Eq. 
(1771) reads (here jc = a%T R /2x) 



P L +xP L + Sn 2 [P L (l+f(t))-f(t)] = 0. 



(83) 



Again we neglect with respect to xPi and P^, obtaining a first-order equation whose solution is 
(in the large x limit) 



Pdt) 



X 



, 2^71; / /8^¥ 2V/27T 
1 + <*> \/ + 



This displays a (quantum) Zeno effect, since for xt 3> 1 one gets 



P L (t) ~ c 



(84) 



(85) 



[compare with (1791) 1 . 

Once again there is a scaling law and one can define a characteristic relaxation time (in natural 
units) 

Tiei!K~xT R . (86) 



maximum of /. 
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Observe that this scaling is at variance with ( T77T ). 

Simulations 
Method 

We will now study in detail the features and results of the integration of the kinetic equation by 
means of a Monte Carlo method already used in the past to study the kinetics of two-level systems 
in nonequilibrium gases.— >^ 

Let us recall the main features of the simulation. Some details have already been given in Sec. 
„ We set Q. = 935 s _1 , a = a R = a L ~ 0.2-r0.4, t~ ! < 2 ■ 10 5 s _1 and N L = N R = 40 energy levels 
in each subspace, with energies given by ©, where n s — 1,...,40, C0l = 1-3 ■ 10 10 s _1 , (Or — 
9.7 ■ 10 9 s _1 . The minimum energy difference AE between the levels is of great importance. One 
can check that with the above-mentioned numerical figures AE/h = 2.8 • 10 9 s _1 and the condition 
(1281) is always satisfied^. The populations dynamics is collected as an average over an ensemble of 
5 ■ 10 3 simulated particles. Since the underlying equations are linear, the particles can be serially 
simulated and the precision of the results sharpened by simply increasing the sample size. The 
simulations provide the time variation of all the elements of the 1 -particle reduced density matrix. 
We constantly checked all the level populations p Hs , 1 < n s < 40, s = L,R, but will only discuss 
in the following the temporal behavior of the total population of the left subspace Pl. The initial 
situation, in all the simulations, is 

Pi L = l, all others = 0, (87) 

so that the initial population is concentrated in the 1 1^) state (the ground state of the left subspace). 

3 The determination of SE = ^^i<m L <N L , \<n R <N R \E mL — E„ R \ for generic N L / R , with E„ s given by (0, poses an 
interesting problem of number theory. However, in our case Nl=Nr= 40 and one can numerically check that the value 
AE/h = 2.8 • 10 9 s~' given in the text is stable against perturbation of (Ol.r of a few percent (well above experimental 
uncertainties). 
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Results 

It is interesting to discuss in more detail some features of the relaxation process and compare them 
to the analytical model proposed in Sec. „ We track the temporal evolution of all the populations 
and try to estimate the speed and the extent at which the levels get populated. Two suitable in- 
dicators are the mean /i = \Xl and standard deviation o = Ol, introduced in (1531) and (l54l) . They 
are plotted in |Figure 6| and |Figure 7| and accurately reproduce the analytical results (1561) and (1571) 



(remember that Dt = xt /Tr). The analytical results are not shown in the graphs, for they cannot be 
distinguished from the numerical ones. 

Notice also in both figures the square-root dependence (1581) for large times t ^> Tr/x ~ 3 ■ 
10 2 Tr. It is worth stressing that this also provides a direct proof that boundary effects, related to 
the finiteness of the number of levels, can be safely neglected for the times considered here. 

oc=0.2 



10 1 



10 L 



— a 2 *x" 1 =32*a/27i 
^ a 2 *x 1 =48*Q/27i 
— a 2 *x 1 =56*Q/27i 




10" 



1 0" 2 1 0" 1 

Time (Rabi periods) 



10 c 



Figure 6: Temporal evolution of the mean ji L introduced in Eq. (I53l) . for <Xr = CCt = a = 0.2 and 

jc = 32, 48, 56. 
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oc=0.2 




10" 3 10" 2 10" 1 10 

Time (Rabi periods) 



Figure 7: Temporal evolution of the standard deviation Ol introduced in Eq. (l54l) . for <Xr = <Xl = 
a = 0.2 and x = 32, 48, 56. 
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We now show how the relaxation of the population depends on the collision frequency X 1 , for 



fixed values of the parameter a = (Xl = OLr. Figure 8 shows the temporal evolution of the relative 



population of the left subspace Pt(t) = 'Ln L Pn L (t) (once again, the analytical results cannot be 
distinguished from the numerical ones and are not shown in the graph). We note that this quantity 
will always eventually tend to its equilibrium value P[ = 1/2, according to (|73T) . However, the 
important point is that by increasing the collision frequency from 300 T^ 1 to 800 T% , the system 
tends to remain in the left subspace for a longer time. This is evident in the plot and is a clear 



manifestation of a QZE. We also notice (although this is not clearly visible in Figure 8 due to 
the scale chosen) that there is always a short-time quadratic region, characterized by a "Zeno time" 
Pl(Q) = — Q m full agreement with Eq. (1681) . The features of this short-time region are independent 
of other parameters (such as a and x),— as can be seen in the figure. Finally, we emphasize that 
x = a 2 T jx ranges between 12 and 32 and is therefore always ^> 1, so that the analysis of Sec. „ 
applies. 

A similar Zeno effect is evident when the parameter a is varied, while keeping the collision 
frequency T _1 constant, as displayed in|Figure 9|(once again, we only display the numerical results, 



for the analytical ones cannot be distinguished). Unlike in the preceding case, where the Zeno 
effect was due to increasing collision frequency, now it is due to increasing collision effectiveness: 
a larger a entails more dephasing and decoherence and, in a loose sense, a better "measurement" 
of the quantum state. The parameter x = a 2 T / ' x ranges between 39 and 72 (^> 1) and one observes 
again the presence of a (parameter-independent) short-time region. 

As the analysis of Sees. n - n shows, the dynamics of the system should be ruled by the scaling 
parameter 

oc 2 T R 

X = DT R = -. (88) 



Figure 10 shows how this scaling law is supported by the results of the numerical simulation. The 
plot shows three sets of curves corresponding to three different values of x. In each set, the values 
of a and x were varied as indicated. Some deviations from the scaling law (of order 5%) can 



be observed and are to be ascribed to the influence of the terms neglected in deriving Eq. ([36 
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cc=0.2 




32 



x 1 =800*n/27i 




* a=0.22 
b a=0.24 
+ — a=0.26 
x a=0.28 
h a=0.30 



0,65 







1 2 3 

Time (Rabi periods) 



Figure 9: Temporal evolution of Pi as a function of a = (Xl = (Xr. For all calculations we set 
x- 1 = 800Q/27Z\ 
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Incidentally, notice again the short-time quadratic behavior. 




0,5 1 1,5 2 

Time (Rabi periods) 

Figure 10: Test of the scaling law (f88l) . Temporal evolution of Pi for different values of x = a 2 T / x: 
three simulations were done with oc^ = = a = 0.2, 0.25, 0.4, respectively, x ranges between 
32 to 56. 

Concluding remarks 

We have studied a Zeno effect in a multilevel molecule made up of 40+40 levels, one of which 
(the ground state of the left subspace) is initially populated, and the evolution towards the right 
subspace is slowed down both because the collisions remove population density "upwards" from 
the left ground state (a classically intuitive process) and because they "dephase" (or analogously, 
make energetically less favorable) the transitions towards the right subspace. The latter process 
is classically less intuitive, but is readily understood if one thinks in terms of quantum transition 
amplitudes (or of the Fermi "golden" rule for a bona fide unstable system). 

It is worth stressing that the general ideas and techniques introduced in this article are valid 
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for any multilevel molecule and any possible level distribution: we focused on the case © only 
for concreteness. Those situations in which (1281) is not valid are very particular cases and their 
analysis, although of interest, goes beyond the scope of this article. 

On the other hand, it is also necessary to emphasize that we neglected temperature effects 
and rapid structural rearrangement phenomena leading to a Boltzmann distribution of the level 
populations. This is a conceptually interesting problem, that involves delicate issues: a sensible 
estimate of the timescales involved in these thermalization processes is a challenging problem that 
requires further investigation. 

We conclude by noticing that the Hamiltonian CD)-© is also relevant for the study of quantum 
chaos and Anderson localization. 35 The analysis of Poissonianly distributed "kicks" © would 
introduce a novel element of discussion in such a context. 
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Appendix 

It is interesting to look explicitly at the derivation of Eq. (1361) from Eq. (1271) . The physical mecha- 
nism at work is the effective decoupling between the fast and the slow modes in (l27l) . Let us start 
from the equation for pi L 2 L , that explicitly reads [here 6>2 L i L = {E% L —E\ L ) /K] 

dpi, 2, „ OCl s 

— = lO>2 L i L pi L 2 L - l£lp\ R 2 L - I— [P2 L 2 L ~ Pl L l L - Pl L 3 L ) 

- 27 (PMl + Pi L 2 L - 2p2 L i L - 2p 2L 3 L + 2pi i2i + Pl L A L ) • (89) 
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When condition (|28l) is satisfied, the first term in the right-hand side dominates over the others and 
one obtains 

^L~ ia>2LlL p lL2u (90) 
which yields a very fast dynamics for the term p\ L 2 L : 

p\ L i L {t) =pi L 2 t (0)exp(iafe £ i i f)- (91) 

The equations for the other off-diagonal components of p are similar. These equations yield very 
rapidly oscillating solutions. 

On the other hand, the dynamics of the populations pi s i s , with s = L,R, and of the coherent 
terms pi L i R is governed by the equations 

dpi, i, n/ , a 2 , a^. a} , , 

-&(P1rIl -PIlIr) + -t(P2z.2 l -PliU) " i —(p2 L l L -pl L 2 L ) ~ ^(p3 L l L + Pl L 3 L ), 



3«, 



P ^ lR = &{P1r1l ~ Pltlfi) + -jr(p2 R 2R - PIrIr) - i-^{p2 R \ R - P\ R 2 R ) ~ ^(p3 R l R + Pl R 

d Phi R jG(n n \ <*L n a l n ,«L n { a R n ,oc L a R 

dt - lil yPl L l L - PIr1r)-^PIlIr-^P1r1l- 1 — P2 l Ir~ 1 — Plz2 R H ~P2 L 2 R 

-^P3l1r-^PU3r- (92) 

It is apparent that no "diagonal" fast frequency ft) is present and these matrix elements evolve over 
timescales % and Q. 1 which are much larger than ft) . Therefore the contribution of all the off- 
diagonal fast terms of the type (|9~T| ) is averaged to zero over the long timescales x and Q -1 , the 
dynamics of the slow and fast terms completely decouples and we get 

~ -^(p lsli -p lilR )+^(p 2L2i -p lilL ), 

^Piffle \ a ff / \ 

— ^ ~ ^(Pl«l L -plz.l«) + Y (p2 « 2 « _Plsls) ' 

^ - zS2 (P1l1l-P1r1rJ- 2tPil1«- 2^Pi«1l- 

(93) 
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Analogously, the evolution equations of the populations p nis = p ms m s read (m s ^ Il,r) 



which are in the form (|29l)-(l36l). Notice the absence of fast and oscillating terms. 
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